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Abstract 

Roy's equations, which incorporate crossing symmetry of the vrvr scattering 
amphtudes, are used to resolve the present ambiguity between two solutions for 
the scalar-isoscalar phase shifts below 1 GeV. It is shown that the "down-flat" 
solution satisfies well Roy's equations and consequently crossing symmetry while 
the other solution called "up-flat" does not and thus should be eliminated. 

1 Introduction 

A good knowledge of the pion-pion scattering is important in studies of different pro- 
cesses in nuclear and particle physics [1]. Construction of phenomenological vrvr ampli- 
tudes requires not only experimental input but also use of theoretical constraints such 
as unitarity, analyticity and crossing symmetry. Direct study of the tttt collisions is 
beyond present experimental possibilities. Phenomenological phase shifts are obtained 
through partial wave analyses of final states in which pions are produced. These anal- 
yses are often model dependent and can sometimes lead to ambiguous results. 

In 1997 a study of the 7r~p| — >■ tt^tt'ti reaction on a polarized target has been 
performed in the 171^^^^ effective mass between 600 and 1600 MeV giving four solutions 
for the TTTT scalar-isoscalar phase shifts below 1 GeV [2]. Using the unitarity constraint 
two "steep" solutions were rejected and the two remaining ones, called "down-fiat" and 
"up-flat" , passed this test [3] . Elimination of this remaining ambiguity is necessary for 
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extension of studies done near the tttt threshold and above using chiral perturbation 
models based on QCD (see for example [4, 5]). 

A steady increase of the tttt scalar-isoscalar phase shifts below 1 GeV can be inter- 
preted as due to the presence of a broad cr-meson [6, 7]. Previous ehmination of "steep" 
solutions excludes a narrow a with a width comparable to that of the p- meson [3]. In 
this channel, near the KK threshold and above it, other resonances exist which could 
be mixed with glueball states predicted by QCD. 

In order to eliminate the above mentioned "up-down" ambiguity one can check 
if the corresponding amplitudes satisfy crossing symmetry. Roy's equations [8] can 
serve as a tool to perform this check. They also correlate different experimental phase 
shifts determined near the tttt threshold and at higher energy, in particular for the 
scalar-isoscalar, £ = 0, J = 0, the scalar-isotensor, £ = 0, 1 = 2, and the vector- 
isovector, i = 1, / = 1, tttt partial waves. These equations have been already applied by 
Pennington and Protopopescu [9] to study different tttt phase shift solutions obtained 
in the early seventies. They were able to resolve a "stccp-fiat" type ambiguity [3] 
present in the phenomenological tttt amplitudes considered at that time. Later two 
high-statistics experiments were performed [10, 11] and, as previously said, the analyses 
of these experiments lead to two plausible solutions [2, 3]. Recently a comprehensive 
analysis of Roy's equations for the tttt interaction has appeared [12]. There, a special 
emphasis is put on the m,r7i range from the tttt threshold to 0.8 GeV. In the present 
paper we pay a particular attention to the range between 0.8 and 1 GeV. In this range 
the largest differences between the "up-flat" and "down-flat" solutions occur, reaching 
values up to 45" (see Fig. 4 in [3] and Fig. 2 here). 



2 Roy's equations 

Assuming analyticity one can write twice subtracted fixed-t dispersion relations [8] for 
the unitary tttt amplitudes. 



T'{s,t,u) = ^Ci/ B^{t) + {s-u)D^{t) 



I'=0 



+ - / ^ f ^ l'' + -^CZ'] Im T^'is', t, u') 

TT J S'^ \S' — S S' —U I 

4^2 \ / 



(1) 



In this equation s, t and u are the usual Mandelstam variables satisfying s+t+u = 4^^, 
where n is the pion mass ; Cj/ and C^l are the isospin matrix elements and 1^^ is the 
unit matrix. The subtraction functions B^' (t) and D^'{t) can be related to the isospin 
and 2 5'-wave scattering lengths Oq and Oq, respectively. 
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Projection of (1) onto partial waves leads to Roy's equations, 



Re //(s) = aQSioSio + alSnSm 

s — 1 1 

12//^ 6 2 ^2) 

O 1 Smax ^ '' 



+ 



r=ov=o ^2 



Here the functions K^l'(s,s') are the kernels and dl{s,Smax) are the so-called driving 
terms. Detailed expressions of these functions can be found for instance in [12, 13]. 
The driving terms contain the low-energy, s' < Smax, contributions from the partial 
waves i' > 2 and the high-energy, s' > Smax, contributions from all the partial waves. 
If one assumes Mandclstam analyticity, the range of validity of Roy's equations (2) 
extends to 68/^2 = (1.15 GeV)^ [8]. 

The partial waves amplitudes //(s) are related to the tttt phase shifts and in- 
elasticities Tj^: 

Each set of experimentally determined phase shifts and inelasticities can serve as input 
to calculate the real and imaginary parts of the partial wave amplitudes. After a 
suitable parameterization, the imaginary parts can be inserted into Roy's equations (2) 
from which one obtains the real parts as output to be compared with the corresponding 
real part input. The quantitative agreement between output and input will be used to 
verify how well a given set of phase shifts satisfies crossing symmetry in a particular 
range of mTTTr where Roy's equations can be applied. 



3 Input amplitudes 

The "down-flat" and "up-flat" data [2] together with the results on the KK phase 
shifts [14] were analysed by us using a separable potential model of three coupled 
scalar-isoscalar channels (vrvr, KK and an effective An system) [15, 16]. This model 
yields particularly good fits to the "down-fiat" solution from 600 to 1600 MeV. The 
"up-flat" data are reasonably well described only above 970 MeV. From 800 to 970 MeV 
the model values are too low in comparison with the "up-flat" data (see Fig. lb of [15]). 
In our calculations we need a faithful representation of both "up-flat" and "down-flat" 
solutions and therefore below 970 MeV we use the following Pade approximation for 
the Sq phase shifts: 

= ni.(F/a.. - 1) 

where k = ^^/V^^TfjiF is the pion momentum and aj{j = 1, . . . , 7, 9) are constant 
parameters. This choice is dictated by the analytical properties of the Jost functions. 
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These constants will be obtained from the best fits to data and to Roy's equations for 
the "up-flat" and "down-flat" separately. In these flts we also use the near threshold 
phase-shift differences, 5q — Si, recently extracted from the high statistics K^^^ decay 
experiment [17]. Using the power expansion of the 5"- wave amplitudes near the tttt 
threshold, 

Ref^{s)^ai + bie + ..., (5) 

we can relate the parameters af] and to the constants aj. One has: Oq = —ai/i 
and bg — —ai/i (0.5/i~^ -|- + a^^ + a^^ — af) — a^fi. The parameters 0:7 and ctg 



are chosen to match the values of 6q following from the three-channel fits at 969 and 
970 MeV [15]. Above 970 MeV we use the fit A to describe the "down-fiat" solution 
and the fit C for the "up-flat" one. These model amphtudes are used up to the Smax 
value even if it exceeds 1600 MeV (see (2)). 

The isotcnsor wave is parameterized within the separable potential model of [7, 15]. 
We use a rank-two potential: 



V^=\p,p') = ^XMp)9^ip'), (6) 

i=l 

where 

are form factors with range parameters p and p' are the pion center of mass momenta 
in the initial and final states, respectively. The threshold parameters Oq and 6q can be 
related to the strength parameters Aj = Aj/ {2j3f). The / = 2 vrvr phase shifts from [18], 
obtained with their method A, serve as input in our fitting procedure. Here we assume 
the isotensor wave to be elastic (?7o = 1) from the tttt threshold to Smax- 

For the P-wave, from threshold to 970 MeV, we use an extended Schenk parame- 
terization as defined in [12] 



tan 



Sl{s) ^Jl-^e{A + Be + Ck^ + Dk^) i tt_IjL j . (8) 




The parameter A is equal to the P-wave scattering length a\ and the parameter B is 
related to the slope parameter b\ = B + AA/{sp — 4:fj,^). The parameter Sp is equal 
to the p-mass squared. Above 970 MeV the P-wave amphtude is represented by the 
K-matrix parameterization of Hyams et al. [19]. The parameters C and D of (8) 
are chosen to match both parameterizations at 969 and 970 MeV. We have checked 
that the near threshold scalar-isoscalar phase shifts determined from the differences 
— ^1 in [17] are insensitive to different P-wave parameterizations. Differences between 
parameterizations of [12] and [19] are smaller than the errors of 5q — 5\. 

The driving terms d[[s, Smax) in (2) are calculated including the contributions 
of 72(1 270) and p3(1690) and the Regge contributions from the Pomeron, p- and 
/-exchanges. We use the Breit-Wigner parameterization of /2(1270) and p3(1690) 
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as described in [2] with masses, widths and tttt branching ratios taken from [20]. The 
range parameters are chosen to be 5.3 GeV~^ and 6.4 GeV~^ for the /2(1270) and 
P3(1690) resonances, respectively (see Table 4 of [19]). The Regge contributions are 
parameterized as in [12] without inclusion of the small w-crossed terms. The Smax limit 
is set to (2 GeV)^ and our results are fairly close to those of [12]. The most important 
contribution in the i = I = channel comes from the /2(1270) resonance. It roughly 
agrees with the result of Basdevant, Froggatt and Petersen [13] where only the /2(1270) 
contribution was considered and s^ax was set equal to (1.46 GeV)^. The dl{s, Smax) I'e- 
sults are much smaller than those of [13] due to the difference in Smax and to the lack of 
the p3(1690) contribution in [13]. In the P-wave there is a strong cancellation between 
the contributions of /2(1270) and p3(1690) leading to very small values, again much 
smaller than in [13]. We have considered the different parameterizations of /2(1270) 
and p3(1690) resonances used by [12]. The changes from the Breit-Wigner form, that 
we used, do not affect the phase shifts in all three partial waves by more than one 
degree below 970 MeV. In the £ — 0, 1 = case the Regge contributions are of the 
order of a few percent of the resonance contributions. For the isospin 1 and 2 they are 
of the same order as the resonance contributions but the corresponding overall driving 
terms are by an order of magnitude smaller than the isospin term. Alternative Regge 
contributions were considered following papers [9], [21] and [22]. These do not change 
significantly the driving terms, the changes being much smaller than one degree in the 
effective mass range up to 1 GcV. 

Our thirteen free parameters, six for the isoscalar ^-wave, four for the isotensor one 
and three for the P-wave are determined through a least square fit to the data together 
with the minimization of the squares of the differences between the input ( "in") and 
output ( "out") of Roy's equations for the three waves. We define 

Xtotal= E [Xlxpi^) + Xloyi^)] ^ (9) 
7=0,1,2 

where 



2 

i=l 



sm 



jsj) - ^j, jsj)) 



2 

(10) 



In (10) (si) and Aipj (sj) represent the experimental phase shifts and their errors, 
respectively. The of the fit to Roy's equations is defined as 



^2 / 7-\ _ fin i^j) ^® font i^j 

XRoyV ) — 2^ 



A/ 



fill 



where sj = [4j + 0.001]/^^ for j = 1, 11 and Su = 46.001yU^. We take a A/ value of 
0.5 X 10~^ to obtain reasonable values of XRoy corresponding to an accuracy of one-half 
percent. Simultaneously we require that the xlxp should not be larger than about 18 
for the fit to the 18 data points between 600 and 970 MeV. We use the CERN MINUIT 
program which provides errors of the fitted parameters. 
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For the isospin wave the number of data points iVo is 24 in (10). It consists of the 
18 "down-fiat" or 18 "up-flat" data points [2] between 600 and 970 MeV and of the 6 
points below 400 MeV from the Kg4 decay experiment [17] . In the isotensor wave we use 
the 12 data points [18] from 350 to 1450 MeV. For the isospin 1 we generate 8 pseudo- 
data points using the X-matrix parameterization between 600 and 970 MeV [19] and 
we choose A(f\ — 2". 

4 Results 

The Fade approximants for the isoscalar 5'-wave amplitude supplemented by the model 
amplitudes [15] for m^rTr > 970 MeV together with the P-wave parameterization and 
that of the isotensor S'-wave as described above, constitute our input to Roy's equations. 

A good global fit to data and Roy's equations can be achieved only for the "down- 
flat" data. The resulting £ = 0, / = phase shifts, plotted as solid lines in Fig. 1 
are compared to the data. The values are summarized in Table 1. For the 18 
points between 600 and 960 MeV the is 16.6 and it is 5.7 for the 6 near threshold 
points. The P®^ degree of freedom on the isoscalar phase shifts is 1.2. Roy's 
equations are very well fulfllled, the differences between the real parts "in" and "out" 
being smaller than 0.8 x 10~^ for the isoscalar wave and smaller than 2 x 10^^ for the 
isotensor amplitude. For the P-wave this difference does not exceed 6 x 10^^. The 
corresponding parameters are given in Table 2. The low-energy parameters Oq, b^, Oq, 
6o and a\, compare well within errors with those of [12] but b\ is smaller by about 
50%. Relative errors of those parameters vary from about 5% for Oq, 6q and a\ to 25% 
for b\. For the other parameters, they are less than 20% with the exception of ^2. 
This parameter is negative, so its influence on the values of 6q is small as can be seen 
from (4) . It furthermore does not appear in the two first coefficients of the low energy 
expansion of tan Sq. One then expects the error of a2 to be large. 

In the "up-fiat" case a good global fit cannot be obtained, as the xlxp '^^ ^^e 18 
data points of [2] between 600 and 960 MeV is as large as 46.4. Including the of 
6.6 from the 6 points of the Ke4 experiment, one obtains the total value equal to 53.0 
(see Table 1). The x^ P^r degree of freedom on the isoscalar data is as large as 2.9 
which shows that this fit for the "up-fiat" case, plotted as the sohd line in Fig. lb, is 
not acceptable. The agreement with Roy's equations is, however, almost as good as 
in the "down-flat" case. In Fig. 2 we plot, together with the data, the solid curves of 
Figs, la and lb representing Roy's flts to the "down-flat" and "up-flat" data. These 
two curves define a band of isoscalar phase shifts fitting well threshold data and Roy's 
equations for two sets of data above 600 MeV. The "down" curve reproduces well the 
"down-flat" data while the "up" curve does not reproduce the "up-flat" data. 

The main differences between the "down-flat" and "up-flat" solutions are for en- 
ergies between 800 and 970 MeV. Below 740 MeV phase shifts of both solutions are 
compatible within their error bars and the "down" curve reproduces well the "up-fiat" 
data (see Fig. 2). We have checked that there is a possibihty to find a good fit to 
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Table 1: values for the different fits; Nj being the number of data and nj that of 
free parameters, A^o — 24, uq — 6, Ni — 8, rii — 3, N2 — 12, 712 = 4. 



"down-fiat" case 
global fit 


"up-flat" 

global fit 


case 

special fit 


I 


1 


2 


1 2 


1 2 




22.3 7.0 
0.1 6.0 


8.1 
0.4 


53.0 8.4 6.9 
0.3 7.4 0.5 


19.0 7.0 8.1 
1.2 10^ 5.5 5.8 


Xtotal 


43.9 




76.5 


1.2 10^ 



Table 2: Parameters from the fit to Roy's equations and to the "down-fiat" data of [2]; 
0:3, 0:7, ccg, B, C, D, Al and A2 are dependent parameters. 



isoscalar S'-wave 


isovector P-wave 


isotensor S'-wave 


Oq 

W 


0.224 ±0.013 
n ocio +0.012 




a\ 
b\ 


(3.96 ±0.24) 10-2 
(2.63 10-3 




2 

4 


I? 


(-3.43 ±0.36) 10" 
(-7.49^1;°^) 10-2 


-2 

/^-^ 


a2 


_4 41 +2-06 
^■^^ -5.05 






30.87 ±0.14 


/^^ 




/3i 


4.88 ±0.25 




0:4 


7.53 ±0.32 












/32 


1.23 ±0.21 




"5 


(2.29 10-2 




















12.0 


















as 


-0.153 




B 


-3.27 10-3 




4 


Al 


0.268 






-0.320 10-=^ 




C 


5.24 10-^ 




6 


A2 


-0.0257 






-0.312 10-^ 


/^-^ 


D 


-2.66 10-^ 
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Roy's equations and to the "up-fiat" data in this hmited range of effective mass. This 
is shown as the dashed fine in Fig. 2. Differences of phase shifts between this fine 
and the "down" one are smaller than one degree up to 840 MeV. Therefore, based 
on the "down-flat" flt, we build up a special Pade flt to the "up-flat" data in order 
to reproduce as well as possible their rriT^T^ dependence. We proceed in the following 
way. Below 600 MeV this flt is constrained to approximate the previously obtained 
"down-flat" isoscalar amplitude and in particular to reproduce its scattering length Cq 
and slope parameter 6q as well as its two values at 500 and 550 MeV. The correspond- 
ing Xea;p(0) ^^.O is good but the XRoy(O) of 1.2 x 10'^ is very large. The differences 
between the "in" and "out" real parts are as large as 0.25 around 900 MeV. So, if one 
tries to improve the flt to the "up-flat" data then one spoils the flt to Roy's equations. 

We have studied the influence of the experimental errors on the vrvr input amplitudes 
by calculating Roy's equations for two extreme isoscalar amplitudes fltted to the data 
points shifted upwards ( "higher-in" ) or downwards ( "lower-in" ) by their errors. Below 
600 MeV these flts were constrained in the same way as in the special Pade flt to the 
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400 600 800 1000 

HI (MeV) 

7T7T ^ ' 




400 600 800 1000 

m (MeV) 

7T7T ^ ' 

Fig. 1: Scalar-isoscalar phase shifts for a) "down- flat" data (full circles) and b) "up- 
flat" data (open circles) [2]. Diamonds denote the Ke4 data [17]. Solid lines represent 
fits to Roy's equations and to data. Dashed fines represent bands of the Pade fits to 
the data of [2]. The dotted hue in b) shows the special Pade fit. 
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400 600 800 1000 

m (MeV) 

7T7T ^ ' 

Fig. 2: Solid lines - fits to Roy's equations and to the "down-flat" data (full circles) and 
to the "up-flat" data (open circles) of [2] as well as to the JCe4 data of [17] represented 
by diamonds. The dashed curve is Roy's flt to the i^e4 data and to the "up-flat" data 
in the restricted range of m-^T, up to 740 MeV. 

"up-flat" data just described above. The corresponding long- and short-dashed fines in 
Figs, la and lb show typical bands delimiting the possible values of the experimental 
data within their errors. The resulting output curves, lower "out" and higher "out", 
of the numerical integrations of Roy's equations for the three partial-wave real parts, 
multiplied by 2A;s~^/^, are displayed in Figs. 3 and 4. In Figs. 3a and 4a we also show 
the input band for the real part of the isoscalar wave, limited by two dashed lines, 
called lower "in" and higher "in" . These dashed lines correspond to the dashed lines 
shown in Figs, la and lb. We do not show the "m" lines in Figs. 3b, 3c, 4b and 4c 
since they are almost indistinguishable from the "oni" lines. 

Below 600 MeV Roy's equations are well satisfied for all waves. At higher energies 
a good agreement is found for the P-wave and 1 — 2 S-wsve. As seen in Fig. 3, in 
the "down-flat" case both "out" curves lie inside the band hmited by the "in" curves 
up to 937 MeV. Above this energy the lower "out" curve starts to lie below the higher 
"in" and at about 950 MeV the higher "out" starts to be above the lower "in" curve. 
The "in" band is then inside the "out" one. It means that there is still a possibility 
to flnd an "out" solution within the "in" band. We can then conclude that within the 
error bars the "down-flat" solution satisfles Roy's equations and consequently crossing 
symmetry. On the contrary, in Fig. 4, for the scalar-isoscalar "up-flat" solution, the 
output band lies outside the input one from 840 to 970 MeV. This eliminates the 
"up-flat" solution as it does not satisfy crossing symmetry in that energy range. 
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Fig. 3: Real parts of the tttt amplitudes (multiplied by 2A;s~^/^) corresponding to the 
"down-flat" data [2] (full circles). In a) input and output bands of the scalar-isoscalar 
real parts are shown. Diamonds denote the K^^ data [17]. Solid and dashed lines in b) 
and c) represent the output bands. In b) triangles denote the isotensor data of [18]. 
Crosses in c) are the isovector pseudo-data calculated from the X-matrix fit of [19]. 
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5 Conclusions 



We have studied the scalar-isoscalar tttt amphtudes constructed from two sets of phe- 
nomenological phase shifts called "up- flat" and "down- flat" [2]. It has been shown 
that only the "down-flat" solution satisfies well crossing symmetry. The conclusion 
about violation of Roy's equations for the "up-flat" data is insensitive to the variation 
of the different parameters in our input: reasonable changes of the i — 1, I — 1 and 
£ — 0, 1 — 2 amplitudes, shift of the upper integration hmit Smax from (2 GeV)^ to 
(1.46 GeV)^, modifications of the /2(1270) and p3(1690) resonance parameters and of 
the Regge amplitudes in the driving terms. All these changes do not lead to complete 
overlap of the input and output bands for the "up-flat" solution. 

Recently a joint analysis of the 5"- wave 7r"'"7r~ data [10, 11] and of the new tt^tt'^ 
data [23], obtained by the E852 Collaboration at 18.3 GeV/c, has been performed [24]. 
It has been shown that, using the one-pion and Oi-exchange model developed in [2], the 
calculated S'-wave intensity of the tt^tc^ system agrees with the measured vr^vr^ intensity 
only for the tttt amplitudes obtained from the "down-flat" phase shifts. So, in [24] the 
"up-flat" solution has also been eliminated and a new "down-flat" set of phase shifts, 
compatible with the tt+tt" and 7r°7r° data, has been found. We have verified that these 
new "down-flat" phase shifts, parameterized as described above, satisfy well Roy's 
equations. 

To conclude, using the theoretical constraints of crossing symmetry, unitarity and 
analyticity, the four-fold ambiguities in the isoscalar 5'-wave tttt phase shifts, in the 
TTTT invariant mass range from 800 to 1000 MeV, have been eliminated in favour of the 
"down-flat" solution. 
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